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A review of measurements of a B is given, representing the status of April 2004. The results prove the energy 
dependence of a s and are in excellent agreement with the expectations of Quantum Chromodynamics, QCD. 
Evolving all results to the rest energy of the Z° boson, the world average of a s (M z o) is determined from mea- 
surements which are based on QCD calculations in complete NNLO perturbation theory, giving 

a s (M z o) = 0.1182 ±0.0027 . 



1. INTRODUCTION 

The coupling constant of the Strong Interac- 
tions, a s , is one of the most fundamental param- 
eters of nature which is to be determined by ex- 
periment. In this review, a summary of the most 
recent measurements of cv s representing the sta- 
tus of April 2004 is given, providing another in- 
cremental update of a more complete and concise 
review pp and of For a detailed introduction 
into the field and for an overview and definition 
of basic concepts, equations and references, the 
reader is referred to 1 1121. 



2. NEW RESULTS 

New or updated measurements of a s are avail- 
able from many classes of high energy particle 
reactions. In the following subsections, the re- 
spective results will be shortly reviewed. 

2.1. Deep Inelastic Scattering (DIS) 

New measurements of a s from inclusive jet 
cross sections in 7p interactions at HERA are 
available from the ZEUS collaboration Pj. Jet 
cross sections and values of a s are presented as 
a function of the jet transverse energy, £^ et , for 
jets with E^ 1 > 17 GeV. The resulting values of 
a s (£^ et ), based on NLO QCD calculations, are 
displayed in Figure 1. They are in good agree- 



ment with the running of a s , as expected by 
QCD, and average to 

a s (M z o) = 0.1224 ±0.0001 (stat.) 

(exp.) 



+0.0022 
-0.0019 
+0.0054 \ 

— 0.0042 t^neo.j, 



(1) 



if evolved to the energy scale of Mzo using the 
QCD f3 function in two-loop approximation, c.f. 
P Equation 7. 

Averaging all measurements of a s from jet pro- 
duction at HERA, which were recently reviewed 
e.g. in 0], results in 

a s (M z o) = 0.120 ±0.002 (exp.) 

±0.004 (theo.), (2) 

unchanged from the previous average 

An update of a s from measurements of polar- 
ized structure functions [5], now containing data 
from the HERMES experiment at the HERA col- 
lider, results in 



a s (M z o) = 0.113 ±0.004 (exp.) 

(theo.), 



-0.009 
-0.006 



(3) 



in NLO QCD. 

A new global analysis using all available preci- 
sion data of deep inelastic and related hard scat- 
tering processes includes recent measurements of 
structure functions from HERA and of the inclu- 
sive jet cross sections at the Tevatron [B]. After 
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Figure 1. a^E^ 1 ) values determined from QCD 
fits of the measured 



dEl 



from ZEUS 0. The 

solid line corresponds to the central value of 
a s (M z o ) determined in [3] , the dashed line repre- 
sents the world average 2 . 



analysis of experimental and theororetica uncer- 
tainties, and restrictions to "safe" fit intervals of 
the input distributions, the authors obtain 



a s (M z o) = 0.1165 



±0.002 
±0.003 



(exp.) 
(theo.) 



(4) 



in NLO QCD. Using NNLO QCD calculations 
wherever available, the same fit gives 



a s (M z a) = 0.1153 ±0.002 (exp.) 

±0.003 (theo.). 



(5) 



The latter result, however, does not relate to com- 
plete NNLO since predictions of jet production 
cross sections and parts of the DIS structure func- 
tions are only available in NLO so far. The previ- 
uosly obtained corresponding result, a s (M z o) = 
0.119 ±0.002 ±0.003 0, was derived without ap- 
plying the new "safe" fit conditions. 

2.2. e + e Annihilation 

A recent summary of a s determinations from 
LEP, at the highest e + e~ collision energies, and 



from previous experiments at the PETRA and 
TRISTAN c + c~ colliders was given in [Hj. Apart 
from the results which were already discussed and 
presented in |1I2| , updates on a s from electroweak 
precision measurements |9I10| . from r lepton de- 
cays and from a recent combination of a s deter- 
minations from hadronic event shape observables 
and jet rates by the LEP QCD Working Group 

The most recent combination of the LEP-I and 
LEP-II electroweak precision measurements of all 
four experiments, in NNLO QCD, resulted in 



a s {M z o) = 0.1226 



±0.0038 (exp.) 

(Mr) 

(QCD) 



+0.0033 
-0.0000 

+0.0028 
-0.0005 



(6) 



from R z = T had /r e = 20.767 ± 0.025, whereby 
the second error accounts for variations of the 
unknown Higgs boson mass between 100 and 900 
GeV/c 2 . The third error comes from a parametri- 
sation of the unknown higher order QCD correc- 
tions, i.e. from variations of the QCD renormal- 
isation scale and renormalistion scheme, see e.g. 

n 

In the same analysis ^U], the fitted leptonic 
pole cross section, cr° = (2.0003 ± 0.0027) pb, 
resulted in 



a s {M z o) = 0.1183 



±0.0030 (exp.) 

(M H ). 



+0.0026 
-0.0000 



(7) 



Since cr9 



has 



p2 

TTTpf and T z ~ T ha d, ere 

lvl z 1 z 

a steeper dependence on a s than has T^^: in 
next-to-leading order, the QCD coefficient C\ for 
T/jad ~ (1 + S„(C„a™)), n = 1,2,3,..., turns to 
2Ci for ct£, C 2 turns to (2C 2 ±Cf ) etc. The exper- 
imental error of a s from ai is thus smaller than 
that from Vhad- However, with increased QCD- 
coefficients Cj, the renormalisation scale uncer- 
tainty also increases, c.f. Equation 13 of 1 , such 
that the QCD uncertainty on a s from ag is ex- 
pected to increase w.r.t. a s from R z . 

A global fit of all LEP data to determine a s 
together with the masses of the Z° boson, of the 
top-quark and of the Higgs boson, gives JOI 



a s (M z o) = 0.1200™i( exp .). 



(8) 
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Figure 2. Summary of measurements of a s {Q 2 ) 
from electroweak precision measurements (filled 
symbols; in NNLO QCD) and from jet and event 
shape observables (open symbols; in resummed 
NLO QCD) in c + e _ annihilation. The dashed 
curves represent the world average of a s (M z o) [5]. 



The latter result is the most precise available from 
combined electroweak fits of the LEP data. 

Finally, using all available data from LEP, from 
SLC and from the Tevatron (i.e. including direct 
measurements of the masses of the top-quark and 
of the W-boson) results in 

a s (M z o) = 0.1186 ± 0.0027 (exp.). (9) 

For the latter two results, there is no additional 
uncertainty due to the unknown Higgs mass. The 
QCD uncertainties for these particular results of 
a s , however, were never determined, and prove 
to be difficult to be guessed due to the unknown 
size of the effective QCD coefficients that enter 
the overall fit. Similar as argued in the case of 
cr», the QCD uncertainty on may be a good 
approximative estimate, but cannot simply be ap- 
plied to other observables - especially if they are 
to be regarded as precision results. 

The combined result of a s from r-decays, in 
NNLO QCD, is [H] 

a s (M T ) = 0.322 ±0.005 (exp.) 

±0.030 (theo.). (10) 



When extrapolated to the energy scale M z o , this 
results in a s (M Z o) = 0.1180 ± 0.0005(exp.) ± 
0.0030(theo.). 

The overall combination of all LEP results on 
hadronic event shapes and jet production rates 
|12| . using resummed NLO QCD predictions, re- 
sults in 

a s (M z o) = 0.1202 ±0.0003 (stat.) 

±0.0049 (syst.). (11) 

This analysis also provides precise values of a s 
in the cm. energy range of LEP-I and LEP-II, 
from 91.2 to 206 GeV. These results, together 
with those from electroweak precision measure- 
ments and from r decays as given above and also 
including a s determinations from lower energy 
e + e _ colliders, are summarised in Figure 2. 

2.3. Hadron Colliders 

Since the last update of this a s review |2j, no 
new measurements of a s from hadron colliders 
were reported. In general, the available results 
from earlier studies, all in NLO QCD, are com- 
patible with but not really competetive to those 
obtained from DIS and from c + c~ annihilation. 
This is due to the sum of systematic uncertain- 
ties, from higher order QCD as well as from the 
available Monte Carlo generators, from underly- 
ing events and beam remnants, and from the en- 
ergy calibration of detectors. Necessary improve- 
ments, like QCD predictions in NNLO and/or in- 
cluding rcsummation, corrections for nonpertur- 
bative effects, improved tools like new and more 
reliable jet algorithms, more data statistics and 
possibly more data at different collision energies 
are, however, underway and should be in place at 
the startup of the Large Hadron Collider. 

3. SUMMARY AND WORLD AVERAGE 

A summary of all significant measurements of 
a s , as discussed in |1I2| and with updates and new 
measurements presented in this review, is given in 
Table 1. 

The values of a s (Q) are presented in Figure 3, 
as a function of the energy scale Q where the 
measurement was carried out. The data provide 
significant evidence for the running of a s , in good 
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Figure 3. Summary of measurements of a s (Q 2 ). Results which are based on fits of a s (M z a) to data 
in ranges of Q, assuming the QCD running of a s , are not shown here but are included in the overall 
summary of a s (-Mzo), see Figure 4 and Table 1. 



agreement with the QCD prediction. 

Therefore it is appropriate to extrapolate all 
results of a s (Q) to a common value of energy, 
which is usually the rest energy of the Z° boson, 
M z o . As described in P , the QCD evolution of a s 
with energy, using the full 4-loop expression 
with 3-loop matching 30 at the pole masses of 
the charm- and the bottom-quark, M c — 1.7 GeV 
and Mb = 4.7 GeV, is applied to all results of 
a s (Q) which were obtained at energy scales Q ^ 
M z o. 

The corresponding values of a s {Mio) are tab- 
ulated in the 4 th column of Table 1; column 5 
and 6 indicate the contributions of the experi- 



mental and the theoretical unceratinties to the 
overall errors assigned to a s (Mzo). All values of 
a s (Mzo) are graphically displayed in Figure 4. 
Within their individual uncertainties, there is 
perfect agreement between all results. This jus- 
tifies to evaluate an overall world average value, 
a^(M z o). As discussed e.g. in pQ, however, the 
combination of all these results to an overall aver- 
age, and even more so for the overall uncertainty 
to be assigned to this average, is not trivial due 
to the supposedly large but unknown correlations 
between invidual results, especially through com- 
mon prejudices and biases within the theoretical 
calculations. 
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Figure 4. Summary of measurements of a s (Mzo). 
Filled symbols represent results based on com- 
plete NNLO QCD calculations. 



For combining all or subsets of the results sum- 
marised in Table 1 into average values of a s (M z o ), 
the same procedures as utilised in pp are being 
used. Averages a^(M z o) for all and for subsets of 
a s -results, together with the corresponding un- 
certainties Acx^ are summarised in Table [21 As 
already discussed in p^, the overall uncertainties 
decrease if the averaging process is restricted to 
those which accomplished a minimum precision, 
i.e. a total error of Aa s < 0.008, while the value 
of a^(M z a) is almost unaffected by such a restric- 
tion - c.f. rows 1 and 2. 

There is a sufficiently large number of redults 
which is based on complete NNLO QCD, such 
that a7(M z o ) can be reliably calculated from this 
subset (see rows 5 to 8 of Tabic |2J . Due to the 
improved completenes of the perturbation series, 



these results are believed to be more reliable and 
better defined than all the others which are com- 
plete to (resummed) NLO. 

The world average of a s (M z o) is finally deter- 
mined from those NNLO results that have total 
errors less than 0.008, namely 



DIS [Bj - SR] 


a B (M z o) = 


n i 01 +0.005 

u.izi_ 009 , 


r decays 


a s (M z o) = 


0.1180 ±0.0030 , 


DIS [v, xF 3 ] 


a s (M z o) = 


n 110+0.007 

U.ll»_Q QQg , 


DIS [e//i,xF 2 ] 


a s (M z o) = 


0.1166 ±0.0022, 


T decays 


a s (M z o) = 


0.118 ±0.006 , 


T(Z -> had) 


a s (M z o) = 


n 1OOfi+0.0058 
u.izzu_ 003g . 



For defining the overall total uncertainty, an "op- 
timized correlation" error is calculated from the 
error covariance matrix, assuming an overall cor- 
relation factor between the total errors of all mea- 
surements. This factor is adjusted such that the 
overall x 2 equals one per degree of freedom. This 
results in the new world average of 



al(M z o) = 0.1182 ± 0.0027 



(12) 



c.f. row 6 of TableEl which is practically identical 
to the final result of the previous summary [2], 
al(M z o) = 0.1183 ± 0.0027. 

The choice of results which contribute to the 
determination of the world average value of 
o7(M z o) is, similar as the calculation and defini- 
tion of its overall uncertainty, arbitrary to a large 
degree. A matter of discussion sometimes is to in- 
clude the overall fit result to the combined LEP 
electroweak precision data (Equation 8) instead 
of, as it is done here, using the result from Rz 
(Equation 6), because the overall fit uses more in- 
formation, the experimental error of a s is smaller, 
and there is no uncertainty on Mh in this case. 
With this replacement, and assuming the same 
QCD uncertainty as for a s from Rz (which may 
not be appropriate, see the discussion in Sec- 
tion 2.2), results in a^(M z o) = 0.1180 ± 0.0029. 
Replacing the result from Equation 6 with the one 
from Equation 9, i.e. with the world fit of elec- 
troweak data, gives aj(M z o) = 0.1177 ± 0.0031. 
In both these cases, 57(M z o) slightly decreases 
and the overall uncertainty increases, but fully 
within the assigned error. The slight increase of 



G 



AcE7 is due to an increased correlation factor of 
0.80 and 0.89, respectively, which is assumed to 
assure that the overall x 2 is unity per degree of 
freedom. For reasons of clarity and of strict def- 
initions of uncertainties, the usage of the result 
from Rz (Equation 6) is preferred. 

The world average of c£(M z o) = 0.1182 ± 

0. 0027 corresponds to the following values of the 
QCD scale A-j^g for different numbers of quark 
flavours Nf, evaluated using the full 4- loop ex- 
pansion of a s and 3-loop matching at the quark 
thresholds (c.f. Equations 4 and 8 of PP): 

A^ 5 = 210±|3 MeV (13) 
A Sf 4 = 294t«MeV. (14) 
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Table 1 



World summary of measurements of a s (status of April 2004): DIS = deep inelastic scattering; GLS- 
SR = Gross-Llewellyn-Smith sum rule; Bj-SR = Bjorken sum rule; (N)NLO = (next-to-)next-to-leading 
order perturbation theory; LGT = lattice gauge theory; resum. = resummed NLO. 





Q 






Aa s (M z ») 






Process 


[GeV] 


a s (Q) 


a s (M z o) 


exp. 


theor. 


Theory 


refs. 


DTS fnnl SFl 


0.7 - 8 




113 + 22i2 

u.110 _ 0.0O8 


±0.004 


+0.009 
-0.006 


NLO 


l5l 


DIS [Bj-SR] 


1.58 


375 + 2-21? 


121 + 

"•J-^-l — 0.009 






NNLO 


rcn 


DIS [GLS-SR] 


1.73 


0.280 + 

xj.^^yj _ 0.O68 


0.112 + 

w — 0.012 


+0.008 
-0.010 


0.005 


NNLO 


HI 


t- dec ays 


1.78 


0.322 ± 0.030 


0.1180 ± 0.0030 


0.0005 


0.0030 


NNLO 


El 


DIS \v\ xF,l 


2.8 - 11 




0.119 + °°nr 
— 0.006 


0.005 


+0.005 
-0.003 


NNLO 


IT51 

1 1 


DIS [e/u; Fol 


1.9 - 15.2 




0.1166 ± 0.0022 


0009 


n 0020 


NNLO 


[Tfil2] 


DIS [e-p — > jets] 


6 - 100 




0.120 ± 0.005 


0.002 


0.004 


NLO 




QQ states 


4.1 


239 + 012 
u.^u^i _ 0.010 


0.121 ± 0.003 


0.000 


0.003 


LGT 




X decays 


4.75 


0.217 ± 0.021 


0.118 ± 0.006 






NNLO 




e+e~ \F2] 

c c L 2 J 


1.4 - 28 




1198 + 2-22^ 

u.±-luu _ 0.0054 


0028 


+ 0.0034 
- 0.0046 


NLO 


12H1 


e^e - [(Tv^hI 

v. L^naQj 


10.52 


0.20 ± 0.06 


130 + 

v.±<jkj _ o.029 


+ 0.021 
- 0.029 


0.002 


NNLO 




e^e - [jets & shps] 


14.0 


0.170 + 2-2?l 
— 0.017 


120 + 2'2n2 
U, - L ^' u — 0.008 


0.002 


+0.009 
-0.008 


resum 


123 
11 


e^e - [jets & shps] 


22.0 


151 + 2-SiS 

U -- L "- L - 0.013 


118 + 2-22! 

1J '- L - LU - 0.008 


0.003 


+0.009 
-0.007 


resum 




e^e - [jets & shps] 


35.0 


145 + °n°rB 
— 0.007 


r, 190 + 0.008 
o.-L^o — 0.006 


0.002 


+0.008 
-0.005 


resum 


123 


e^p — [(T^^^fl 


42.4 


0.144 ± 0.029 


0.126 ± 0.022 


0.022 


0.002 


NNLO 


|23ll| 


e + e _ [jets & shps] 


44.0 


— U.UUo 


n 190 + 0.008 
— u.Uub 


0.003 


+0.007 
-0.005 


resum 


EZJ 

1 1 


e + e~ [jets & shps] 


58.0 


0.132 ±0.008 


0.123 ±0.007 


0.003 


0.007 


resum 


m 


pp -> bbX 


20.0 


n U r + 0.018 
U.140 _ 0.019 


0.113 ±0.011 


+ 0.007 
- 0.006 


+ 0.008 
- 0.009 


NLO 


m 


pp, pp -> 7X 


24.3 


f) 10c + 0.012 
U.lOO _ Q.008 


n llf) + 0.008 
U.11U _ 0.005 


0.004 


+ 0.007 
- 0.003 


NLO 


m 


o-(pp -> jets) 


40 - 250 




0.118 ±0.012 


+ 0.008 
- 0.010 


+ 0.009 
- 0.008 


NLO 


m 


c+c- [r(Z -> had)] 


91.2 


n 199A+ 0.0058 
U.1ZZU_ 0.0038 


fl 199fi+ 0.0058 
U.1ZZU_ 0.0038 


±0.0038 


+0.0043 
-0.0005 


NNLO 


ma 


e + e~ seal. viol. 


14 - 91.2 




0.125 ±0.011 


+ 0.006 
- 0.007 


0.009 


NLO 




e + e~ 4-jet rate 


91.2 


0.1170 ± 0.0026 


0.1170 ±0.0026 


0.0001 


0.0026 


NLO 




e + e~ [jets & shps] 


91.2 


0.121 ±0.006 


0.121 ±0.006 


0.001 


0.006 


resum 


□ 


e + e~ [jets & shps] 


133 


0.113 ±0.008 


0.120 ±0.007 


0.003 


0.006 


resum 




e + e~ [jets & shps] 


161 


0.109 ±0.007 


0.118 ±0.008 


0.005 


0.006 


resum 


□ 


e + e~ [jets & shps] 


172 


0.104 ±0.007 


0.114 ±0.008 


0.005 


0.006 


resum 


n 


e + e~ [jets & shps] 


183 


0.109 ±0.005 


0.121 ±0.006 


0.002 


0.005 


resum 


n 


e + e~ [jets & shps] 


189 


0.109 ±0.004 


0.121 ±0.005 


0.001 


0.005 


resum 


n 


e + e~ [jets & shps] 


195 


0.109 ±0.005 


0.122 ±0.006 


0.001 


0.006 


resum 


m 


e + e~ [jets & shps] 


201 


0.110 ±0.005 


0.124 ±0.006 


0.002 


0.006 


resum 


m 


e + e~ [jets & shps] 


206 


0.110 ±0.005 


0.124 ±0.006 


0.001 


0.006 


resum 


m 
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Table 2 

Average values of a^(M z o) and averaged uncertainties, for several subsamples of the available data. The 
result printed in bold-face is taken as the new world average value of 57(Af z o). "— F 2 " denotes omission 
of the result from F2 structure function fits JHI ■ 











opt. corr. 


overall 


uncorrel. 


row 


sample (entries) 


a;(M z o) 




correl. 




1 


all 


(32) 


0.1190 


0.0038 


0.69 


0.0009 


2 


" Aa s < 0.008 


(23) 


0.1192 


0.0034 


0.63 


0.0010 


3 


all - F 2 


(31) 


0.1196 


0.0043 


0.71 


0.0010 


4 


" Aa s < 0.008 


(22) 


0.1198 


0.0039 


0.68 


0.0010 


5 


NNLO only 


(9) 


0.1182 


0.0031 


0.68 


0.0015 


6 


" Aa s < 0.008 


(6) 


0.1182 


0.0027 


0.61 


0.0015 


7 


NNLO - F 2 


(8) 


0.1195 


0.0042 


0.75 


0.0019 


8 


" Aa s < 0.008 


(5) 


0.1196 


0.0038 


0.74 


0.0020 


9 


NLO only 


(23) 


0.1197 


0.0044 


0.69 


0.0011 


10 


" Att s < 0.008 


(17) 


0.1199 


0.0040 


0.65 


0.0012 


11 


DIS only 


(6) 


0.1169 


0.0033 


0.76 


0.0018 


12 


DIS - F 2 


(5) 


0.1185 


0.0060 


0.77 


0.0031 


13 


e + e~ only 


(22) 


0.1199 


0.0044 


0.77 


0.0012 



